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Abstract: Within the framework of a (l+l)-dimensional model which mimics evolution and 
scattering in QCD at high energy, we study the influence of the running of the coupling on 
the high-energy dynamics with Pomeron loops. We find that the particle number fluctuations 
are strongly suppressed by the running of the coupling, by at least one order of magnitude 
as compared to the case of a fixed coupling, for all the rapidities that we have investigated, 
up to y = 200. This reflects the slowing down of the evolution by running coupling effects, in 
particular, the large rapidity evolution which is required for the formation of the saturation front 
via diffusion. We conclude that, for all energies of interest, processes like deep inelastic scattering 
or forward particle production can be reliably studied within the framework of a mean-field 
approximation (like the Balitsky-Kovchegov equation) which includes running coupling effects. 
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1. Introduction 

In the recent years, one has assisted at two major directions of progress in our understanding 
of the QCD dynamics at high energy, both aiming at improving over the previous status of the 
theory, as encoded in the BaHtsky-Kovchegov [1,2] (BK) and Jahhan-Marian-Iancu-McLerran- 
Weigert-Leonidov-Kovner (JIMWLK) equations [3-5]. One of these directions started with the 
recognition [6-8] of the fundamental role played by gluon-number fluctuations (or 'Pomeron 
loops') in the QCD evolution with increasing energy, which led to more complete evolution 
equations (the 'Pomeron loop equations') [8-14] and to a surprising link to problems in statistical 
physics [7], which is rich in consequences [8,15-19]. The other line of research concentrated on 
the inclusion of next-to-leading order (NLO) effects in the non— linear BK equation [20-27], and 
very recently culminated in an explicit calculation of the running coupling effects associated with 
fermion loops within perturbative QCD [28-33]. Both types of effects — Pomeron loops and 
NLO corrections — turn out to be effects of 0(1), which modify in a dramatic way our previous 
expectations based on the BK-JIMWLK equations. Yet, so far, these two types of effects have 
never been considered together in a unified framework (in particular, their mutual influence has 
never been addressed), mainly because of the complexity of the Pomeron loop equations, which 
are quite difficult to deal with already at leading-order. 

In this paper, we shall for the first time study the effects of the running of the coupling 
on the high-energy evolution with Pomeron loops, within a one-dimensional model proposed 
in Ref. [34] which captures in a simplified form the relevant dynamics in QCD. For this model, 
that we shall extend here to the case of a running coupling (in such a way to mimic the one-loop 
running coupling of QCD), we shall provide a rather detailed numerical analysis of the 'dilute- 
dense' scattering process similar to deep inelastic scattering (DIS) in QCD. Our results point 
out towards a very interesting, and rather unexpected, conclusion: The 'Pomeron loop' effects 
(i.e., the effects of particle-number fluctuations) are strongly suppressed by the running of the 
coupling — by at least one order of magnitude as compared to the case of a fixed coupling — , 
in such a way that they remain negligible up to the highest energy that we have investigated, 
which corresponds to a rapidity Y = 200. This finding has a corollary of great practical interest: 
it implies that, for studies of the high-energy dynamics (say, in the energy range at LHC) one 
can reliably resort on appropriate mean-field approximations, so like the BK equation properly 
generalized to include NLO (or, at least, running coupling) effects. 

Our numerical analysis is further corroborated by analytic estimates which help clarifying 
the mechanism responsible for the suppression of fluctuations: this is not the fact that, because 
of the running, the value of the coupling is effectively smaller, as one might naively think; rather, 
this is related to the fact that the running of the coupling drastically slows down the evolution. 
With a running coupling, not only the saturation momentum grows much slower with the energy, 
but the same is also true for the "BFKL diffusion"' [35] , which is the mechanism through which 
the saturation front (say, for the gluon occupation number, or for the dipole scattering amplitude 
in DIS) evolves towards its asymptotic, 'geometric scaling', shape at high energy. Namely, the 
diffusive radius grows with Y like Y^^^ (rather than the usual y^/^ behavior at flxed coupling), 
and such a rise is too slow for the front to reach its asymptotic shape within the energy range 
under consideration. Rather, the front preserves a pre-asymptotic shape, which is not favorable 
for the growth of fluctuations. Hence, during this whole 'pre-asymptotic' evolution (that we 
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found to extend up to F = 200 at least), the DIS dynamics is the same as predicted by the 
corresponding mean-field approximation (with running coupling, of course). 

Although obtained in a specific model, we are confident that these conclusions should apply 
to QCD as well, for several reasons: First, as just mentioned, these results are supported by 
analytic estimates which are essentially identical in QCD and in the model under consideration. 
Second, this model is truly similar to QCD: it has been constructed [34,36] by requiring boost- 
invariance together with an evolution law inspired by the gluon evolution in the context of 
JIMWLK equation (see Sect. ^ below for details). While simpler than the original JIMWLK 
equation, because of the absence of color degrees of freedom and the reduction to only one spatial 
dimension (which, physically, plays the role of the gluon transverse momentum), the model is 
at the same time more general, in that it allows for particle-number fluctuations. In particular, 
the mean-field approximation to this model turns out to be quite similar, even quantitatively^, 
to BK equation of QCD. Third, this model belongs to the universality class of the statistical 
process known as 'reaction-diffusion' [37,38], so like the high-energy QCD dynamics itself [7]. 
Although, strictly speaking, this universality has been so far established only for the case of 
a fixed coupling, we expect all the processes in this class to respond in a similar way to the 
inclusion of a running coupling. 

When applied to QCD, our conclusions could in particular explain why phenomenological 
analyses inspired by the mean-field (BK) dynamics were relatively successful in describing the 
small-x experimental situation at HERA, despite the fact that the proton wavefunction was 
expected to develop strong correlations via Pomeron loops in the course of the evolution : such 
correlations, which would transcend any mean-field description, are in fact suppressed by the 
running of the coupling. Most significantly, the phenomenon of geometric scaling — a hallmark 
of the BK dynamics [21-24,39] which seems to be weU verified by the HERA data [40-42], 
but which would be washed out by the fixed-coupling evolution with fluctuations [7, 8, 17] — 
is in fact resuscitated by our present analysis. But this analysis also shows that, precisely due 
to running coupling effects, the window for strict geometric scaling is drastically reduced as 
compared to the respective mean-fleld estimate at fixed coupling [21,22]. (The width of this 
window is controlled by the 'diffusive radius' alluded to above; see Sect. ^ for details.) This 
conclusion is consistent with previous studies of the BK equation with running coupling [22-24] , 
and also with the phenomenological analyses of the HERA data [43-45] , which need to include a 
substantial amount of scaling violation in order to accurately describe the experimental results. 

This paper is organized as follows: In Section 2, we briefly review the construction of the 
one-dimensional model [34] and show how this can be accommodated with the running of the 
coupling. Section 3 offers a summary of known results concerning the dynamics with Pomeron 
loops at fixed coupling, as well as the mean-field dynamics with running coupling. Based on 
such previous results, we formulate our theoretical expectations for the full problem (Pomeron 
loops and running coupling) at the end of that section. These expectations are then confronted 
to numerical results in Section 4, which is devoted to an extensive analysis of the model. Section 
5 contains our conclusions and some perspectives. 



^In the sense of generating similar numerical values for the anomalous dimension, the saturation exponent, 
etc.; see Sect. H. 
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2. Model description 



The construction of the model from the underlying physical assumptions — boost invariance, 
multiple scattering in the eikonal approximation, and evolution law a la JIMWLK (meaning 
that the probability for emitting a new particle in one step of the evolution depends in a non- 
linear way upon the density of particles created in the previous steps, and saturates at high 
density) — has been described in detail in a previous publication [34] (see also Refs. [36, 46] 
for earlier version of the model, in zero spatial dimensions), so here we shall only describe the 
generalization of this model to the case of a running coupling. 

The "hadronic systems" which undergo evolution and scattering are two systems of classical 
particles ('the right and left mover') distributed in one spatial dimension ('the x-axis') which is 
transverse to the collision axis. Prom the point of view of QCD, the position x of the particles 
along this transverse axis corresponds to the logarithm of the transverse momentum of a gluon 
(or, within the QCD dipole picture [47], to the logarithm of the inverse size of a dipole). The 
total rapidity gap between the incoming systems is equal to Y, and we choose the 'laboratory 
frame' in such a way that the right (left) mover carries a rapidity Y — Yq (respectively, — io)- 
The scattering is assumed to be elastic, and the corresponding 5-matrix {S)y to be real (as 
typically the case in QCD at high energy). Needless to say, {S)y must be independent upon the 
rapidity divider Yq, i.e., upon the choice of the frame. 

The particle configurations in the two systems are created via high-energy evolution (from 
their respective rest-frames to the laboratory frame), which is a stochastic process. Accordingly, 
these configurations are themselves random, and can be described in the language of probabil- 
ities: we shall denote by P^[n{x^),Y — Yq] and, respectively, Pi,[m(xi^),YQ] the probability 
densities to find given configurations in the two systems. Note that a configuration is described 
as a function n{x) (the density of particles at point x), and the probabilities introduced above 
are functionals of this density, as well as functions of the rapidity variable. 

The average ^-matrix is then computed as the following, functional, integral: 

{S)y = JvnVmPK[n{xK),Y -Yo]PL[m{xL),Yo]S[n{xK),m{xL)], (2.1) 

which represents the average over all possible configurations of the 'event-by-event' S'-matrix 
S[n(xYi),m{xi^)] associated with a given pair of configurations. In turn, the latter is given by^ 



S[n, m] = exp 



dxRda;Ln(xR)m(xL) In ^(xrIxl) 



(2.2) 



where o"(xr,|xl) = 1 — T{xYi\xi^) is the 5-matrix for the scattering of two elementary particles 
of logarithmic sizes xr, and xl, and t(xr|xl) the corresponding T-matrix (a real quantity, in 
between and 1). 

Consider now the evolution of any of the incoming systems with increasing rapidity. An 
'evolution step' corresponds to a small increment dY, and we shall assume that only one extra 
particle can be emitted in such a step. Moreover, the evolution law will be taken to be quite 
simple: when one particle is emitted, the final configuration consists in the same particles as 



^Perhaps it is easier to understand this expression if we start from a discretized version where one has S[n, m] 



Ylij cr^ ' ' and then take the continuous hmit to arrive at Eq. (2.2) 
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the initial configuration plus an additional particle of arbitrary size. We can quantify this 
assumption in terms of the following master equation (with 6xz = S{x — z)) : 



Mn{x)-dx,]P[n{x)-6x,,Y] - / /,[n(x)] P[n(x), F]. (2.3) 



dP[n{x),Y] 

dy ~~ 

The quantity fz[n{x)] is a "deposit" rate density, that is fz[n{x)]dzdY is equal to the probability 
that we will find an extra particle with logarithmic size in the interval (z, z + dz) after one 
evolution step dY, given that the initial configuration was n{x). The interpretation of the two 
terms in Eq. ( |2.3| ) as gain (for the positive one) and loss (for the negative one) terms is then 
straightforward. 

Lorentz invariance (i.e., the condition that Eq. ( |2.1] ) be independent of Yq) then requires [34] 
the deposit rate to be proportional to Tu][n{x)] — the T-matrix for the scattering of a particle 
of logarithmic size to off a system at a given configuration n{x) : 



Tu,[n{x)] = 1 - exp 
We fix the proportionality constant by choosing 



dxn{x) lncj(tt;|x) 



(2.4) 



fM^)] = ^^- (2.5) 

where a{z) is the running coupling in the problem. This precise dependence of the deposit rate 
upon a is chosen in such a way to recover the correct powers of a in the evolution equations to 
follow, in agreement with perturbative QCD. It is important to notice that Eq. ( |2.5| ) together 
with Eq. ( |2.4| ) for the scattering amplitude imply that the deposit rate density is in general a 
highly non-linear function of the particle density. In turn, this means that the extra particle at 
z is emitted coherently from all the preexisting particles in the system. 

To completely specify the model, we also need to specify the form of the elementary particle- 
particle scattering amplitude T{x\y) and the x-dependence of the running coupling a{x). The 
former will be chosen as 

T{x\y) = a{x)a{y) exp(-|a; - y\) = a{x)a{y)K{x, y), (2.6) 

which mimics the corresponding quantity in QCD, i.e. the elementary dipole-dipole scattering 
amplitude at zero relative impact parameter. Indeed, to lowest order in perturbative QCD 
(where the coupling is fixed), the dipole-dipole amplitude (at zero impact parameter) is ap- 
proximately proportional to a^r^/r?,, where r< = min(ri,r2) and r> = max(ri,r2); this is in 
agreement with Eq. ( |2.6| ) once one recalls that x corresponds to the logarithm of the inverse 
dipole size. Of course, after including running coupling effects, which are part of the NLO cor- 
rections, the actual QCD amplitude becomes more complicated than its toy model counterpart 
in Eq. (|2.6| ). Here however we shall stick to Eq. ( |2.6D , which is general enough to illustrate our 
point, and has also the advantage to be symmetric under the interchange x <^ y. 

As far as the running coupling a{x) is concerned, we would naturally define it as 1/ (3x, where 
/3 is a fixed number — the analog of the one-loop beta-function in QCD. Such a definition would 
introduce a divergence at x — > 0, similar to the infrared problem of perturbative QCD. However, 
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in the context of the non-hnear evolution of interest for us here, such a divergence is in fact 
innocuous, except perhaps at low energies, because the saturation momentum introduces an 
effective infrared cutoff, which grows rapidly with Y (see Sect. ^). Still, in view of the numerical 
simulations, it is necessary to regulate the singularity of the running coupling at x — > 0; the 
onset of saturation will then guarantee that the physical results are independent upon the 
precise regularization prescription (at least, at sufficiently large Y). In what follows, we choose 
to "freeze" the coupling when entering the "non-perturbative" region at x < 0. (Note that our 
present conventions are such that x = corresponds to a transverse momentum fc^ = Aq^-jj-, in 
QCD.) Namely, we shall require 



1/ (5x for X » 1 

0(1) for x~0(l) (2.7) 
Qo < 1 for X = — oo, 



with q(x) a monotonically decreasing function of x. A precise implementation of this running 
will be presented in Sect. |^. 



Let us now return to Eq. (2.5) for the deposit rate and study its behavior in the limits of 



low and high particle density. We easily find that 

\ f dx a{x)K{z,x)n{x) when n(x) <C l/a(x)a(z) for all x 
[1/0(2;) when n{x) ^ 1 / a{x)a{z) for some x. 

Therefore, when the system is dilute the extra particle at z is (approximately) incoherently 
emitted from any of the preexisting particles. This limit, where the deposit rate is proportional 
to the particle density, and of 0(as), is analogous to the QCD dipole picture [47]. When the 
particle density starts to increase the extra particle is emitted coherently, and in the limit where 
the system becomes very dense the deposit rate saturates, at a value of 0{l/as)-, and thus 
becomes independent of n. This is again analogous to the situation encountered in QCD, where 
the JIMWLK equation leads to the saturation of the gluon emission rate [48]. 

To conclude the presentation of the model, let us show how to construct evolution equations 
for observables. In what follows, we shall consider two types of observables: the particle density 
and the scattering amplitude for a collision in which the left mover consists in a fixed number 
of particles (meaning that the whole evolution up to rapidity Y is given to the right moving 
system^). These observables can be cast in the form 

{0)y = j VnP[n{x),Y]0[n{x)], (2.9) 

with the understanding that P[n(x), Y] = Pji[n{x),Y] for the scattering problem. Differentiating 
this equation w.r.t. Y and making use of the master equation ( |2.3| ) one finds 

d{0)Y _ 
dY 



dz (/.[n(x)] {0[n{x) + 5,,] - 0[n{x)]}). (2.10) 



^ Given the boost-invariance of our formalism, this choice brings no loss in generality in so far as the evolution 
is concerned, but merely amounts to specifying the initial conditions for the left mover. 
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In particular, the evolution of the average particle density is governed by (dropping the index 
Y and using subscripts instead of arguments from now on) 

with the natural interpretation that the rate of change in the average particle density at x is 
equal to the average deposit rate density at x. In general, this is a complicated non- linear 
equation, which involves (upon the expansion of the r.h.s.) all the fc-body particle densities, like 



the particle pair density (fc = 2), etc. But at low density, Eq. ( 2.11 ) reduces to a linear equation 

d(n^. 



dY 



Ja,K^,{n,) for n « l/a^ , (2.12) 



analogous to the running coupling version of the BFKL equation [35] for the dipole density [47]. 

Consider now the T-matrix = 1 — Sx for the scattering of a single left moving particle of 
logarithmic size x off a generic right moving system. For a given configuration of the latter, this 
is given by the general expression Eq. (|2.2| ) with m{xi^) = 5(xl — x). Making use of the general 



evolution equation ( 2.10 ) we easily arrive at 



dY 



a. 



j Kxz{T,{l-Tx)) (2.13) 



an equation analogous to the first Balitsky equation (extended to running coupling). Once again, 
this is not a closed equation, as it involves the T-matrix (T^T^) for the scattering of a projectile 
made with two particles. It thus becomes necessary to write down the corresponding evolution 
equation, which is similarly found as 

^^^^ = a,, j Kxz{T,Ty(l - Tx)) +ay j Ky,{T,Tx{l - Ty)) 

z z 

+ a^ayj a,KxzKy,{T,{l - T^)(l - Ty)), (2.14) 

2 

This equation is analogous to the second Pomeron loop (PL) equation [8,9]. Notice that it is 
the last term, proportional to a^, in the r.h.s. of the above equation which distinguishes this 
equation from the corresponding one in the analogous B-JIMWLK hierarchy. At a first glance, 
this last term seems to be suppressed with respect to the first two terms when counting powers 
of a. However this is not true, since this term becomes a dominant one for T ~ a^. Let us briefly 
discuss two limits of the hierarchy starting with the above equations. While this hierarchy, as 
it stands, is not consistent with a solution of the factorized form {TT...T) = {T){T)...{T), it 
becomes so when we keep only the terms which are explicitly proportional to a. This is the 
mean field approximation (MFA) in which the whole hierarchy reduces to the factorized form of 



the first equation ( 2.13 ), that is, 



gy --ax I Kxz{T,){l-{Tx)), (2.15) 
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which is analogous to BK equation with running couphng. If we further drop the nonUnear term 
proportional to (T)^ we arrive at the toy-model analog of the BFKL equation for the dipole 
scattering amplitude. Notice that in this case the running coupling stands outside the integral in 
the r.h.s., to be contrasted with Eq. ( p.l2 ) for the particle density where one needs to integrate 
over the z-dependence of the coupling. As it should be clear from the previous examples, in the 
present model it is the parent particle which sets the scale for the argument of the coupling. 



3. From known results to expectations 



In this section, we shall summarize the known results concerning the fixed coupling case (with or 
without fiuctuations) as well as the running coupling case in the mean-field approximation. We 
shall specialize these results to the one-dimensional model under consideration, but one should 
keep in mind that similar results hold in QCD as well, under the corresponding approximations. 
Based on such previous results, we shall then formulate some theoretical expectations for the 
full problem — the stochastic evolution with running coupling. 

The most important property of the non-linear evolution is the emergence of the saturation 
'momentum' XgiY), which is the scale separating between the dense and dilute regions, and 
which increases with Y. This scale is an intrinsic property of an evolving particle system, but 
it also acts as the unitarization scale for the scattering between that (evolving) system and an 
external dipole. For sufficiently large Y, the scattering amplitude T(x, Y) takes the form of a 
front which interpolates between a strong scattering (saturated target) region at x < Xg, where 
T = 1, and a weak scattering (dilute target) region at x > Xg, where the amplitude decreases 
exponentially. (For the stochastic evolution, this front picture holds event-by-event.) One can 
then define the saturation line Xs{Y) as the 'position of the front', that is, the line along which 
the amplitude is constant and of 0(1), e.g. T{x = Xs{Y),Y) = 1/2. 

(i) In the fixed coupling case, the energy-dependence of Xg is rather well understood, both 
in the MFA, and in the full evolution including fluctuations. 

(i.a) In the mean-field case, and for sufficiently large Y, the 'front velocity' 



>^s{Y) 



1 dxsiY) xi7s) 



a 



dY 



7s 



2-/saY 



3^3 
2aY 



(3.1) 



slowly approaches the constant, asymptotic, value Aq = x{ls)hs = 3\/3 from below. Here 
xil) = 2/(1 — 7^) is the eigenvalue function of the linear equation ( |2.12 ) for fixed coupling and 
7s = l/\/3 solves the equation x{ls) ^ IsX'ils) = 0. Note that, for fixed coupling (FC), it is the 
combination aY which is the natural evolution 'time'. Furthermore, the amplitude in the front 
region (the region ahead of the saturation line but relatively close to it) is obtained as 



r(x, Y) = ci{x - Xs + C2) exp 



-7<;(X - Xs) 



[X 



2x'^aY 



(3.2) 



valid for 1 <^ x — Xg ^ 2x'saY. In Eq. ( |3.2| ), ci and C2 are unknown constants of 0(1), and 
Xs = x"{ls) = 27. We further notice that within the more restricted window 



1 < X 



< x^miY) = ^2x'laY , 



(3.3) 
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Figure 1: The (average) saturation scale and the corresponding velocity in the evolution at fixed cou- 
pling, as obtained via the numerical study of the one-dimensional model, for 3 values of a. The results 
of the full evolution (thin lines) are compared to the respective predictions of the MFA (thick lines). 



where the diffusion term in the exponent in Eq. (^]^) can be neglected, the amplitude exhibits 
geometric scaling, i.e. it depends upon the kinematical variables x and Y only via the difference 
z = X — Xs{Y). 

For even larger x, such that x — Xs ^ 2x'sO(Y , the amplitude shows 'color transparency', i.e. 
it exhibits a faster, exponential, decrease with x, which is the same as the long-range decay of 
the elementary amplitude in Eq. (|2.6|) : T(x) oc e^^. 

(i.b) After 'mcludvag fluctuations, one must distinguish between the event-by-event front, 
corresponding to an individual realization of the stochastic evolution, and the statistical ensemble 
of fronts, which determines the average quantities. 

Up to fluctuations, the asymptotic velocity is the same for all fronts, and thus is also 
the same as the average (asymptotic) velocity. But in contrast to the mean-field scenario, 
this asymptotic value is now approached much faster (exponentially in aY), and its value is 
substantially smaller than Aq- Numerical calculations show that the stronger a is, the larger is 
the reduction in the value of A^ (see Fig. |l|; our procedure for extracting the front position from 
the numerical results will be explained in Sect. An analytic estimate for this reduction is 
known only in the limit where a is extremely small, in which case one finds 

Ao 2x(7s)ln2a2 2ln^ ' ^ ' ' 

This reduction in the value of the speed is in turn associated with a profound change in the 
nature of the front, which is now compact. This property is most transparent in the discretized 
version of the model, as appropriate for numerical simulations: then, the number of particles 
per bin is discrete, and thus it cannot be at the same time non-zero and less than one. Hence, 
the (event -by-event) front associated with the particle occupation number n(x, Y) cannot have 
an infinitely extending exponential tail, but rather it suddenly ends at some 'foremost occupied 



- 8 - 



bin' xq. The front for T{x,Y), on the other hand, is not truly compact, since an exponential 
tail T{x) ~ 0^6"^^"^"^ is generated at x > xq by the scattering off the particles located around 
Xq. But in this case the 'compactness' refers to the fact that the front width xq — Xg cannot 
be larger than L ~ (1/7^) log(l/a^), which is the distance over which the amplitude decreases 
from its value T = 1 at saturation down to a value T ~ 0{a'^). When this maximal value is 
reached, one says that 'the front has been formed'. Prior to that, the width of the front grows 
via diffusion, VY, so the front formation requires a typical 'formation time' 

L2 In^a^ ^ ^ 

To summarize, for fixed coupling and Y > Iform; the event-by-event front for the scattering 
amplitude has the following, approximate, shape, up to fluctuations : 



T{x,Y) = < 



1 for z < 

A z e"^^^ for l<z<L (3.6) 
5a2e-(^-^) for z > L, 



with constant factors A and B of 0(1). The precise interpolations between the shown regimes, 
as well as the shape of the front at earlier stages, Y < Ifomn when the front has not yet fully 
formed, are not under analytic control. Note also that, under the present approximations, the 
event-by-event amplitude shows geometric scaling: T{x,Y) ~ T{x — Xs{Y)). 

It turns out that the shape of the individual fronts is also important for the statistical 
properties of the ensemble of fronts which is obtained by repeating the same evolution a large 
number of times. The fronts which compose this ensemble differ from each other via their 
respective front position Xs, which now is a random variable. The correspondence [7] with the 
reaction-diffusion process in statistical physics suggests that, to a very good approximation, the 
distribution of Xg at y is a Gaussian, with an expectation value (xg) and a dispersion cr^ which 
both rise linearly with Y [15,49]: 

{xs)iY) = XgaY , a\Y) = {xD - {xgf = DaY . (3.7) 

An analytic estimate for the front diffusion coefficient D is available only for a — > 0, in which 
case one finds [15] a similar kind of logarithmic behavior as in Eq. ( |3.4| ): 



with constant do. The very slow, logarithmic, convergence of the results in Eqs. ( |3.4| ) and ( p.^ ) 
to their respective mean-field limits (k = 1 and D = 0) as a — > reflects the strong sensitivity 
of the fixed-coupling evolution to fluctuations. This is conflrmed by numerical simulations (in 
particular, for the model under consideration [34]), which also show that, for more interesting 
values of the coupling a (say a = 0.1 0.5) the coefficients D and 1 — k = (Aq — As)/Ao are 
sizeable numbers, of 0(1) (see Figs. |l]and|2|). Hence, as manifest in Fig. |2|, the dispersion cr'^{Y) 
rises quite fast with Y, which in turn has profound consequences on the shape of the average 
amplitude {Tx) [7, 8] : when averaging over all the fronts in the ensemble, the geometric scaling 
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Figure 2: The front dispersion in the fixed-couphng evolution for 3 vahies of a; the respective values of 
the diffusion coefficient D can be read off the figure on the right. Note the 'formation time' Iform (which 
increases when decreasing a) during which the dispersion remains negligible. 



property characteristic of the individual fronts, cf. Eq. (p.6|), is washed out by dispersion and 
eventually replaced, when a'^iY) ^ 1, by a new form of scaling known as 'diffusive scaling' [17] : 
namely, {Tx)y scales as a function of the single variable {x — Xs{Y))/^/Y. 

Moreover, the numerical simulations demonstrate a delay in the onset of dispersion (in the 
plots in Fig. |2|, this delay is especially visible in the curves corresponding to smaller values 
of a), which can be roughly identified with the formation time Iform- This relation between 
the formation time for the individual fronts and the onset of dispersion is crucial in order to 
appreciate our forthcoming results at running coupling, and can be explained as follows: The 
growth in a"^ is associated with the rare fluctuations in which a new particle is produced relatively 
far away ahead of the tip of the front (see, e.g., the discussion in Ref. [15]). Such a fluctuation 
is most effective when the splitting takes place near the tip of the front: first, this enhances the 
probability that the daughter particle be produced further away ahead^; second, the occupation 
numbers near the tip are low, of 0(1), so the fluctuations are relatively important there. 

However, this whole picture starts to apply only after the front has been formed, i.e., 
when the tail of the front n{x,Y) has developed the characteristic 'anomalous dimension' 7^ 
all the way down to its tip. Assume that, at Yq = 0, one starts with the initial condition 
n(x,0) = no0(xo — x), with no = 0(1)- Then, in the originally occupied bins at 2; < xq, 
the occupation numbers will rise very fast with Y, as n{x,Y) ~ nge^"^ [34], until they reach 
a large, 'saturation', value, of 0{l/a); this happens around Kjat — (l/2a) ln(l/a), which at 
weak coupling is parametrically smaller than Yform- Thus, during the intermediate evolution at 
Ysat ^Y^ YfoT-ra, the front n{x, Y) shows a very abrupt fall off (much faster than the exponential 
law Q-'ys(^-^s) eventually reached at y ~ yform) from a saturation value rigat ~ l/a down to 
n = 1. Such an abrupt fall off is not favorable for developing fluctuations, since the front has no 

''it has been argued in Ref. [15] that this probability behaves like P{S) oc e'"'^"^ where S is the distance from 
the tip of the front Xtip « + L to the place where the new particle is generated. 
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distant tip with occupation numbers of 0(1)); hence, during this interval in Y the front evolution 
remains 'mean-field'-hke. Fluctuations start to significantly develop only for Y > Iform- 

(ii) We now turn to the running coupling case, which so far has been investigated only 
at the mean-field level, i.e., at the level of the BK equation in QCD, and at different levels of 
sophistication for implementing the running coupling effects [20-32] . For the present model, the 
corresponding, non-linear, equation in shown in Eq. ( 2.15| ). 



Such previous studies have shown that, for the purpose of computing the high-energy asymp- 
totic of the saturation scale, one can assume that the front speed is determined locally by the 
corresponding fixed coupling result [21,50]; that is, 



lf?.A,„(.,).^^..(r).^?^r, (3.9) 



where we have used the asymptotic FC result in Eq. ( |3.lD together with a{xs) = l/(/3xs). This 
simple argument is justified since the relevant dynamics takes place within a limited distance 
ahead of Xs, of the order of the 'diffusive radius' Xdis- The latter rises with Y, but this rise is 
considerably slower than that of the saturation scale itself; thus, the relative width Xdis/xg of 
the active region decreases with Y. In fact, the rise of the diffusive radius with increasing energy 
is much slower with RC than with FC [22,23,26] (see below), and this difference turns out to 
have crucial consequences for the physics of fluctuations, as we shall discover later on. 

Based upon the above considerations, it is possible to perform a more accurate study [22-24] 
of the solution to the mean-field equation ( |2.15p , in which the argument x of the running coupling 
is self-consistently expanded around Xs{Y). One thus finds (up to a multiplicative constant) 

(3.10) 



T(x, Y) ~ expl-Ux - X,)] r'/'Ai L + ^^^^) eJip 
where the saturation line Xs(Y) is more accurately determined by 



3x"(7.)t 



, dxs Xofl \Ci\Ds 1 , 



Here we have defined r = (2Ao / (3){Y + Yo), with Yq an unknown constant. Furthermore, 
Ds = {x"{'ys)/['^xi7s)]Y^'^ ~ 1-65, Ai is the Airy function, and ~ —2.34 is the location 
of its rightmost zero. As compared to the earlier estimate in Eq. ( |3.9| ), the expression ( 3.11J ) 



also provides the first pre-asymptotic correction to the front velocity, to be compared to the 
corresponding FC result in Eq. ( |3.1D . 

The above formulae show that it is natural to interpret the variable r cx ^/Y/P as the 
'evolution time' with RC. For the purposes of the numerical calculations in the next section, it 
is also useful to define a 'velocity' with respect to this natural 'evolution time'; we shall write 

Vs = , ~ V 2Ao 1 : 777^ ■ (3.12) 

d7(yZ5) ^ '\ 4 r2/3; 

At large Y, this velocity approaches the constant value vq = \/2Xo which is independent of /?. 
The expression ( p.lO| ) for the amplitude is valid in the interval 

l<x-x. «^r2/3, (3.13) 
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which defines the 'front region' for RC. When one restricts oneself to the narrower window 

1 < X - X, « Xdiff(l^) = , (3.14) 

then one finds that the amplitude shows a 'geometric scaling' behavior, 

T{x,Y) ~ (x - Xs + c)exp[-'ys{x - Xs)], (3.15) 

formally similar to the FC case. Once again, the upper bound of the geometric scaling region 
defines the diffusive radius. As anticipated, with RC, this radius grows very slowly with Y : 
XdiS ~ T^^^ ~ y^/^, to be compared to the much faster increase at FC : x^is ~ 
Eq. ( p. 3D . Some numerical results for the mean-field evolution with running coupling (RC) 
(which will confirm the above analytic estimates) will be presented in Sect. ^. 

We note that in order to obtain the above scaling in Y of the diffusion radius, it was 
rather crucial that we did not set the scale determining the argument of the coupling to be the 
saturation scale. (Such an approximation would not be correct and would lead to a diffusion 
radius proportional to ^/T ~ yi/^ which is smaller than the actual one.) In Eq. ( ^l0|) , we 



have two sources of diffusion, namely, the Airy function, for which x^iff ~ t^/'^ ~ Y^^^, and 
the Gaussian function, for which Xdig ~ t^/^ ~ Y^^^. In the overall product, it is the smaller 
of the two diffusive radii that controls the physics; indeed, the diffusion radius plays the role 
of the effective phase space for evolution, and therefore the smaller it becomes, the larger the 
preasymptotic corrections are^. Indeed, had we set the coupling equal to l/(3xs, we would have 
found a smaller correction, of order 1/r^, to ( p. 11 ). 



(iii) Based on the previous discussion, we are now prepared to formulate some theoretical 
expectations for the results of the full problem, which includes both running coupling and particle- 
number fluctuations. In the next section, these expectations will be confronted to the actual 
numerical results. 

Let us start with the asymptotic regime, where the situation looks conceptually simpler. 
(We shall later specify when we expect this asymptotic regime to install.) Then, by the same 
argument as discussed in relation with Eq. (3.9), we expect the rate of change in both the average 



saturation scale and the dispersion be governed by the respective results at fixed coupling, cf. 
Eq. ( |3.7| ) evaluated with the local value of the coupling a{{xs))- That is: 

^^-^oa{{xs)) = ^, ^^Daiixs)). (3.16) 

From the FC case, one should recall that the coefficients k and D in the above equations depend 
upon a (the stronger the larger is a), and hence they become y-dependent in this RC case. 
However, for sufficiently large values of Y, the coupling a{{xs)) becomes arbitrarily small, and 
then the estimates in Eqs. ( p.4| ) and (|3.8| ) become reliable. These formula show that, in this 
asymptotic regime, the dependence of the coefficients k and D upon a (and hence upon Y) is 
merely logarithmic and thus can be neglected in integrating Eq. ( p. 161) . One then finds: 



51- 



the FC manifest on Eq. (pAh 



Notice that the relative correction in (3.11) is of order r ' ~ x^ia- This "inverse-square law" also holds in 
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where, strictly speaking, the factors k and D have a weak dependence upon Y, and the inte- 
gration constant Yq is not under control. Moreover, k, should be strictly smaller than one, and 
it should slowly approach one from below when Y ^ oo. In the same limit, D should slowly 
approach zero, and the front velocity Vg defined as in Eq. (3.12) should approach the value 



vq = \/2Ao; this is the same limiting value as in the mean-field case, but the convergence to- 
wards it should be much slower with fluctuations (because of the slow convergence of k towards 
1) than in the MFA. 

While the above considerations look reasonable indeed, they do not tell us how fast is this 
universal high-energy regime approached when increasing Y. This question is even more crucial 
at RC than it was at FC, since with RC, the high-energy limit y — > oo is tantamount to the 
weak coupling limit a ^ 0. Hence, the most interesting phenomena may be concentrated in the 
early (or 'pre-asymptotic') stages of the evolution, where the coupling is stronger. 

In order to answer this question (at a qualitative level, at least), we need to study the front 
'formation time' with RC. This is determined by the same general argument as for FC, namely 
this is the rapidity evolution Yform which is required for the diffusive radius to become as large 
as the width L of an individual front. From Eq. ( 3.15| ) (which remains approximately true for 



the event-by-event fronts even in the presence of fluctuations), L is evaluated as: 



11 I Y 

L ~ -In^— ~ -In-, (3.18) 



and hence it is slowly increasing with Y. By also using the RC-version of the diffusive radius 



Xdiff, as given by the upper limit in Eq. (3.14), one finds 



. L ^ ^ . ^ (A) , (3^,9) 

This is, strictly speaking, a transcendental equation (since L itself depends upon Iform) albeit 
only slowly, cf. Eq. ( 3.18| )), that we shall not attempt to explicitly solve here, since our argument 



is at best qualitative. Rather, for our purposes, it suffices to notice that, with RC, the formation 
time scales like yform ~ PL^, and hence it is parametrically larger than in the FC case (where 



we have seen that iform ~ {^/oi)L , cf. Eq. (3^)). In a previous discussion in this section, we 



have argued that the formation time for the individual fronts also acts as the onset time for 
the growth of dispersion. We thus conclude that, with RC, the onset of the fluctuation effects, 
and thus of the universal, asymptotic, behavior in the high-energy evolution, should be strongly 
delayed. Our previous, analytic, estimates are too crude to more quantitatively characterize 
this delay. We therefore turn to a numerical analysis, with results to be presented in the next 
section. 



4. Numerical results 



In this section, we shall numerically study the evolution of the particle distribution in a 'hadronic' 
system described by our model (with running coupling and fluctuations), and also the scattering 
between this system (the 'target') and a simple, unevolving, projectile, which consists in a single 
particle of variable 'size' x. We shall consider initial conditions where the particle density shows 
a plateau at Iq = : n(x,0) = noQ{xo — x), with no chosen in such a way that the original 
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occupation numbers (the number of particles per bin) be of 0(1). In practice, we divide the x axis 
into bins of width Ax = 1/8 and start the evolution with 2 particles per occupied bin, meaning 
riQ = 16. (We have checked that our results are insensitive to the discretization prescription.) 
This choice introduces a dissymmetry between the target and the projectile already in the 
initial conditions: the target is relatively denser, although still far away from saturation; this 
dissymmetry will be, of course, further amplified by the subsequent evolution. 

In our simulations, we shall pay special attention to ensure that the final results are insen- 
sitive to the prescription used to 'freeze' the coupling at negative values of x (cf. Eq. ( |2.7|) ). 
This is necessary in order to unambiguously distinguish between the results of a perturbative 
evolution with running coupling (which is the problem of interest for us here) and the 'non- 
perturbative' phenomena associated with the freezing of a. Clearly, this is an issue only for 
the evolution in the early stages : for sufficiently large values of Y, the saturation scale is large 
too, Xs{Y) ^> 1, and acts as a 'hard infrared cut-off' which removes any sensitivity to the 'non- 
perturbative' region^ at x < 0. However, if the freezing is important in the early stages, then 
this early evolution is 'fixed-coupling'-like, and thus it is relatively fast. Hence, the results of 
the evolution accumulated in these early stages can be numerically important and difficult to 
disentangle from those of the running coupling evolution at later stages. To avoid this, we shall 
make sure that, in the whole range of x covered by the evolution, the running coupling takes its 
perturbation shape a(x) = l//?x to a very good accuracy. 

In practice, we shall use the following prescription for the running coupling, which provides 



a smooth interpolation (actually, a family of such interpolations) to the behavior in Eq. (2.7): 



Ci(x) = 7 -, --; , (4.1) 

^ ' /?cln(e^'/c + ei/"o/3c) ^ ^ 

where < c < 1. For x > l/(/?ao), this running coupling becomes independent of its value in 
the infrared, oq. The additional parameter c allows us to control how fast is this transition, 
from a frozen coupling to a perturbative one: the smaller c, the sharper the transition. We have 
checked that, with c = 0.1, there is no influence of freezing as soon as x ^ 2/{f3ao). 

In the numerical calculations to be presented below, we have chosen uq = 0.7, c = 0.1, and 
a value for xq (the upper bound of the initial plateau at y = 0) which is related to (3 in such a 
way that the condition xq > 2/(/?ao) be satisfied. We have considered several values for /3 (the 
parameter which controls the speed of the perturbative running), but payed special attention to 
the case /5 = 0.72, which is similar to the one-loop beta function of QCD (with Nf = 3 flavors 
of quarks). We have performed systematic calculations up to a maximal rapidity Yma,x = 200. 

The fundamental quantity in our analysis is the position of the front (or 'saturation scale') 
Xs{Y), which can be defined as a line of constant amplitude for the dipole scattering — say, 
T(x = Xs{Y),Y) = 0.01. It turns out, however, that a more convenient definition in practice is 
the following one'' (T(x, Y) is the event-by-event amplitude) : 

Xs{Y) = xs{0) + J [T{x, Y) - T(x, 0)] , (4.2) 



^This should be contrasted to the Unear, BFKL-hke, evolution, where the sensitivity to the non-perturbative 
region at low momenta persists up to arbitrarily large Y, because of the 'infrared diffusion'. 

^We have checked that our numerical results are unchanged if one instead defines using T{x — Xs) — est.. 
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Figure 3: The front dispersion in the evolution with running coupUng, for Y < 200 and 3 values of (3. 



which exploits the fact that for a given event the amplitude satisfies geometrical scaling, i.e. it 
is a function only of the combined variable x — Xs{Y), cf. Eq. ( |3.6|) . (This property is not exact, 
but it is satisfied for all values of x that give the dominant contribution to Eq. (|4.2| ).) This 
definition together with Eq. ( [4.2[ ) implies the following representation for the front dispersion 

a'iY) ^ - {x,f = I [{T^Ty) - {T,){Ty)] , (4.3) 

xy 

where the expectation values in the r.h.s. are computed at rapidity Y. 

In practice, of course, we shall extract Xs{Y) for each individual front via Eq. ( |4.2| ), and 
then construct statistical quantities, so like (xg), cr'^, and also the 3rd order cumulant, C3 = 
(x^) — 3{x'^){xs) + 2{xs)^, via a statistical analysis over the ensemble of events. (Notice that C3 
would exactly vanish if the event distribution was Gaussian.) 

We are now in a position to present our numerical results. The main conclusion emerging 
from these results can be succinctly formulated as follows: the 'Pomeron-loop' effects (i.e., the 
influence of the particle-number fluctuations on the evolution towards saturation) are strongly 
suppressed in the evolution with running coupling, and remain negligible for all the rapidities 
of interest. The physical origin of this suppression is the slowing down of the evolution by the 
running of the coupling and, more precisely, the large 'formation time' required for the formation 
of the front and the onset of fluctuations. This conclusion, and its above physical interpretation, 
are supported by the ensemble of the numerical evidence we now present. 

(i) The strong suppression in the dispersion o"^ and the 3rd cumulant C3 

In Fig. 1^ we have displayed our numerical results for the dispersion a'^{Y), for rapidities 
Y < 200 and two values of the parameter f3. The dispersion grows with Y, as expected; 
moreover, as visible in the right hand plot, this rise is roughly consistent (at least for sufficiently 
large values ofY^ 100) with the y/Y-law expected according to Eq. ( |3.17| ). 
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Figure 4: Left: A comparison between FC and RC results for the dispersion. Right: The ft'ont diffusion 
coefficient with running couphng, as extracted from a fit to the numerical results. 



What is perhaps less expected, and thus surprising at a first fight, is the extremely small 
magnitude of the measured dispersion, as compared to the respective fixed coupling results 
(compare in this respect with Fig. §). One might think that this reduction in cj^ is due to the fact 
that, with RC, the coupling is effectively weaker, but this is actually not true. To demonstrate 
this, we have compared in the left hand side of Fig. Q the dispersion cr^ produced in the RC 
run with (5 = 0.72 and, respectively, the FC run with a = 0.1. This particular value a = 0.1 
is appropriate for our argument since it is close to the smallest value of the running coupling 
reached in the corresponding RC simulation; namely, this is the same as a{{xs)) = \/l3{xs) 
with (xs) measured at y = 75 (see below). Yet, as manifest on Fig. ^ the total dispersion 
accumulated in the RC evolution is tremendously smaller (by almost two orders of magnitude) 
than the corresponding result with FC. 

A similar situation occurs for the 3rd order cumulant C3, as illustrated in Fig. |5[ We thus 
conclude that the strong suppression of the effects of fluctuations that we observe in our results 
cannot be imputed to the fact that the coupling is effectively weaker in the RC scenario than in 
the FC one (for a same value of a at y = 0). 

Incidentally, since the y-dependence of our results for cr^(y) appear to be rather well 
described, at least at large Y , by the expression in Eq. (3.17), it is interesting to perform a fit 



based on this expression and thus extract the value of the front diffusion coefficient from the 
numerical results. We have repeated this fit for several values of /3, for rapidities within the range 
200 < Y/ (3 < 300, with the results shown in the right hand side of Fig. As indicated there, 
these fits suggest the scaling law D{(3) oc 1//3, for which we have no fundamental understanding. 

(ii) The average front position and velocity 

Further evidence for the suppression of fluctuations comes from a study of those quantities 
which exist already at mean field level, like the (average) position of the front (xg) and the 
corresponding velocity. Our respective numerical results are displayed in Fig. |6|, for both the 
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Figure 5: The 3rd cumulant C3 in the evolution with running couphng, for Y < 100 and 3 values of /?. 
For comparison, the fixed coupling result corresponding to a = 0.1 is also shown; note that, in order to 
fit inside this plot, the FC result has been divided by a factor of 20. 




Figure 6: The average front position (left) and the corresponding velocity (right) for the running 
coupling evolution: the results of the full evolution, including fluctuations (thin lines) are compared to 
the respective mean-field results (thick lines). 



complete evolution and its mean field approximation. These results should be compared to the 
corresponding ones at FC, in Fig. |l[ In the latter, the difference from the mean field behavior 
is manifest, and is increasing with Y. In Fig. ^, on the other hand, there are only minor 
differences between the full evolution and the MFA: the respective curves fall almost on top of 
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each other, for all values of Y . (The small visible differences are in fact consistent with the 
magnitude of the dispersion reported in Fig. ^.) In particular, on the right hand figure, one 
sees that all the velocity curves remain well below the limiting value vq = \/2Xq ~ 3.22 which 
is expected at asymptotically large Y both in the MFA, cf. Eq. ( |3.12| ) , and in the full evolution 
with fluctuations. These results have the following implications: 

(a) For the considered range in rapidities, both types of evolution (with or without fluctu- 
ations) are still far away from the asymptotic regime. 

(b) In the full evolution including fluctuations, the rapidities covered by our simulations 
are too small for the fronts to be fully formed. (Since, if the fronts were fully formed at some 
intermediate value of Y, then their subsequent evolution would have been quite different from 
the mean field case.) In particular, the fact that the measured front velocity remains signifi- 
cantly smaller than the asymptotic value vq is not to be attributed to 'cut-off' effects due to 
fluctuations, but rather to preasymptotic effects, which for Y < Yform are the same with or 
without fluctuations. 

Conclusion (b) above can be reformulated by saying that the front 'formation time' Iform 
should be larger than, or marginally comparable to, the maximal rapidity which is covered 
by our analysis. This is also consistent with the previous results concerning the smallness of 
fluctuations: as explained in Sect. |3|, the fluctuations cannot significantly develop until the 
individual fronts have been fully formed. 

The fact that the front 'formation time' Iform can be relatively large at RC should not be 
a surprise, in view of the discussion in Sect. ^. What is quite surprising, however, is that this 
quantity can be that large, namely of O(IOO) or even larger. (In fact, in our simulations, we have 
never seen any sign that the fronts have been formed, until the highest rapidities.) To further 
consolidate this discovery, we performed another numerical test, which directly measures the 
front region, and thus the formation time: 

(iii) The reduced front: event~hy~event geometric scaling 

By inspection of the analytic results in Sect. H (see, especially, Eqs. (^^) and (3.10) there). 



one sees that a study of the reduced front Ti.ed(a^, Y) = e"'''^^~^''^T{x, Y) can give us access to the 
width of the geometric scaling window, for both fixed and running coupling. Indeed, for x ^ Xs, 
which includes the saturation region and the transition region around x ~ x^, where the shape 
of the amplitude is not analytically known, Tj-ed{x,Y) vanishes exponentially, whereas for very 
large x, such that x — Xg ^ a^jifr, it rapidly vanishes once again, due to diffusion. On the other 
hand, within the window for geometric scaling, Tred oc {x — Xs + c) grows linearly with x — Xs, 
cf. Eq. ( p. 15 ); hence, the width of the geometric scaling window can be identified as the region 



of linear increase for T,.e(^{x,Y). 

We have measured the reduced fronts with both EC and RC, with the results shown in Fig. 
1^. For fixed coupling (figure on the left), there is a clear difference between the MFA and the 
complete evolution including fiuctuations: whereas in the former case, the scaling window keeps 
growing with Y, and quite fast (recall that the respective diffusive radius is expected to grow 
like VY, cf. Eq. (|3.3| )), in the presence of fluctuations, the growth of the scaling region is similar 
in the early stages, but then it saturates, around Y ~ 20, at a maximal value which defines the 
width of the front. The fiuctuations in the tail of the front at large x — Xg are clearly visible. 
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Figure 7: The reduced front as a function x — Xg for different values of F, with of without fluctuations (in 
the fluctuation case, we show the respective average quantity). Left: FC evolution with a = 0.1 (from left 
to right: r = 5, 10, . . . 50). Right: RC evolution with (3 = 0.72 (from left to right: Y = 10, 20, . . . 100). 
Notice the different scales used in the two figures. 



Moving to the RC case (right figure), the situations looks very different: the curves cor- 
responding to mean-field and, respectively, fluctuations are now very close to each other, no 
fluctuations are visible in the tail, and, besides, the width of the scaling region is much smaller 
than for FC, as expected (since with RC this region rises only like Y^'^ , cf. Eq. (^) ). The 
fact that, in the fluctuation case, the width of the scaling region keeps growing with Y confirms 
that the front has not been formed, in agreement with the previous results in this section. In 
fact, the maximal extent of the geometric scaling window as seen on this plot is about x — — 5 
for Y = 100, which is rather small — at most, marginally comparable to our previous estimate. 



Eq. (3.18), for the front width L evaluated for a = 0.1. 



(iv) Approximate geometric scaling for the average amplitude 

Since our numerical results show very small dispersion, cf. Fig. |3|, it is quite clear that the 
shape of the average amplitude (T(x))y remains close to that of any of the individual fronts 
which compose the ensemble, and hence to the shape of the mean-field front. We thus expect 
{T{x))y to exhibit approximate geometric scaling, according to the pattern described below 
Eq. ( p.lOj ). In order to check this, we have displayed in Fig. ^ the average amplitude emerging 
from our RC simulations (with /3 = 0.72) as a function of the scaling variable z = x — (xg), for 
various values of Y up to l^ax = 100. From this figure, one sees that the various curves fall 
almost on top of each other for z ^ 5, but they start to deviate from each other for larger values 
of z, although this deviation remains quite small. This behavior is indeed consistent with both 
the analytic considerations in Sect. ^ and the previous numerical results in this present section. 
Namely, the value z ~ 5 for the upper bound of the geometric scaling window is in agreement 
with our former results for the reduced front, as shown in the right plot in Fig. |7[ Also, the fact 
that the dispersion between the various curves (at z > 5) is decreasing when increasing Y, as it 
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can be checked via a close inspection of Fig. ^, is in agreement with our theoretical expectation 
that the violation of geometric scaling proceeds via BFKL diffusion (cf. Eq. ( p. 10 )). 

The quality of the geometric scaling behavior^ visible in Fig. ^ is considerably lower than 
for the corresponding mean-field results at FC, but much better than it would be for the full 
FC evolution, including fluctuations. (For Y = 100, the stochastic FC evolution with a = 0.1 
would generate a dispersion ~ 10, which is large enough for the geometric scaling to be 
completely washed out in the average amplitude, and replaced by diffusive scaling; cf. Fig. |2| 
and the discussion in Ref. [34].) In fact, the amount of geometric scaling violation in Fig. |8| is 
indeed comparable to the one observed in the small-x data at HERA [40-42]. 

We conclude this section with a lesson from our analysis which may shed some light on 
a rather surprising recent finding in Ref. [32]: namely, in that work, an unexpectedly large 
'anomalous dimension' was reported, 7 ~ 0.85, based on an exponential fit to the fronts generated 
by the BK equation with running coupling. This value 7 ~ 0.85 is significantly larger than 
the saturation anomalous dimension 7^ ~ 0.63 which is expected to control the decay of the 
amplitude in the geometric scaling region. Or, as we have seen in both Fig. ^ (right figure) 
and Fig. ^, the window for true scaling is in fact quite narrow: z ^ 5 at Y = 100; outside this 
window, the amplitude decays much faster, because of the diffusion term. It is therefore likely 
that the result reported in Ref. [32] is the consequence of enforcing an exponential fit over a 
relatively large window in z, much larger than the actual window for geometric scaling. 



This quality could be more precisely characterized with the method proposed in Ref. [42]. 
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5. Conclusions and outlook 

In this paper we have presented the first analysis of the consequences of the running of the 
couphng on the high-energy evolution with Pomeron loops, in the context of a simple one- 
dimensional model which captures the relevant dynamics in QCD. We have found that the 
Pomeron loop effects are strongly suppressed by the running of the coupling, up to the high- 
est energies that we have investigated — which go well beyond the energies of interest for the 
phenomenology of QCD. The main reason for this surprising behavior is the fact that, in the 
presence of a running coupling, the diffusive radius grows very slowly with rapidity. There- 
fore, the wave fronts preserve a pre-asymptotic shape which is not favorable for the growth 
of fluctuations. During this pre-asymptotic evolution, which in our simulations extends up to 
Y ~ 200, the dynamics is similar to the respective prediction of the mean-field approximation 
(with running coupling, of course). In particular, due to the suppressed dispersion, approximate 
geometric scaling is preserved for the average scattering amplitude. The window for strict geo- 
metric scaling is rather narrow, but an approximate scaling behavior extends well outside this 
window. We believe that this should suffice to explain the respective scaling observed in the 
HERA data, but a firm conclusion in that sense would deserve a dedicated study. 

In the course of our analysis, we have performed various numerical tests to verify that our 
conclusions are insensitive to the details of the numerical procedure (e.g., the discretization of 
the spatial axis), to the prescription used the 'freeze' the running of the coupling in the infrared, 
and also to the specific implementation of a running coupling in the context of the model (in 
addition to our main choice, Eq. ( |2.6| ), for the running of a"^, we have also performed calculations 
with some other prescriptions, and obtained similar results). 

It is also useful at this stage to summarize the limitations of our present analysis, and 
thus pave the way towards further studies: This study has been restricted to asymmetric initial 
conditions (dilute projectile vs. dense target), where the projectile corresponds to a dipole in 
QCD. This setup covers those physical problems which, in QCD, admit a dipolar factorization — 
chiefly among them, deep inelastic scattering, but also particle production at forward rapidities 
in hadron-hadron (pp, pA) collisions. However, our current model can accommodate arbitrary 
initial conditions, including more symmetric ones, and it would be interesting to study the 
influence of the initial conditions on our main conclusion (the suppression of Pomeron loop 
effects by the running of the coupling). This could shed more light on the physical mechanism 
responsible for this suppression and, in particular, consolidate the argument about the front 
formation time that we have proposed, and which seems to be supported by the present analysis. 

Furthermore, it would be important to clarify the model-dependence of our conclusions, 
if any. If the extension of the present analysis to real QCD (which would require solving the 
Pomeron loop equations of Refs. [8, 10]) looks still prohibitive for the time being, it would 
be nevertheless interesting to apply a similar study to other one-dimensional models of the 
'reaction-diffusion' type, and thus check whether the universality expected for such models still 
persists after the inclusion of running coupling effects. This may have physical consequences for 
other kinds of problems, say in the framework of statistical physics, where what we call 'running 
coupling' in the context of QCD would actually correspond to an inhomogeneous medium, in 
which the rates for particle splitting and merging scale with x like (powers of) 1/x. 
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